Invisible decays of ultra-high energy neutrinos 
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Gamma-ray bursts (GRBs) are expected to provide a source of ultra high energy cosmic rays, 
accompanied with potentially detectable neutrinos at neutrino telescopes. Recently, IceCube has 
set an upper bound on this neutrino flux well below theoretical expectation. We investigate whether 
this mismatch between expectation and observation can be due to neutrino decay. We demosntrate 
the phenomenological consistency and and theoretical plausibility of the neutrino decay hypothesis. 
A potential implication is the observability of majoron-emitting neutrinoless double beta decay. 

PACS numbers: 13.15.+g 13.35.-r 14.60.Pq 14.60.St 95.30. Cq 95.85.Ry 98.70.Sa 



The source of ultra high energy cosmic rays remains a 
mystery. In gamma-ray burst (GRB) models such as the 
fireball model cosmic-ray acceleration should be accom- 
panied by neutrinos produced in the decay of charged 
pions created in interactions between the high-energy 
protons and 7-rays [IJ. Recently the Ice-Cube collab- 
oration reported an upper limit on the flux of energetic 
neutrinos associated with GRBs almost four times be- 
low this prediction [2]. Here we speculate on the plausi- 
bility of the decay hypothesis as a possible explanation 
for the mismatch between observation and expectation. 
The most attractive possibility involves invisible decays, 
which have been considered theoretically since long time 
ago [3]-[7]. These arise in models with spontaneous vi- 
olation of ungauged lepton number [H [8]. A natural 
scenario to test neutrino stability are astrophysical ob- 
jects [9HT2]. In particular, limits on Majoron couplings 
from solar and supernova neutrinos have been obtained 
in Refs. [13 (for non astrophysical constraints for ex- 
ample from OuPP searches, see [ 14HT7] ). Moreover, re- 
cent results from Pierre Auger Observatory (PAG) \1S\ , 
ANTARES [19^ and IceCube [20] have put strong con- 
straints on the neutrino flux coming from distant ultra 
high energy (UHE) neutrinos. Some mechanisms to ex- 
plain a possible deficit in the observed neutrino flux have 
been suggested [21| [22j and the neutrino decay hypothe- 
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sis has been recently revisited [23] and an estimate of the 
relevant GRB neutrino lifetime ranges been given in [24j . 

Here we explore the phenomenological plausibility and 
theoretical consistency of the decay hypothesis within 
a class of low-scale SU(3)c (8) SU(2)l ^ U(1)y seesaw 
schemes with spontaneous family-dependent lepton num- 
ber violation. We show that the required neutrino de- 
cay lifetime range hinted by the non observation of 
UHE muon neutrinos is theoretically achievable for the 
majoron-emitting neutrino decays and, moreover, consis- 
tent with all existing phenomenological constraints. 

The decay rate Ui Uj -\- J in the rest frame of Ui is 
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where Vi and Vj are active neutrinos and J is a massless 
or very light majoron. Taking mj = 0, we can estimate 
the decay length for a relativistic neutrino as given by 
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For typical AGN distances we obtain the required values 
of gij for a given neutrino mass which would cause 
decay before reaching the detector. In Fig. ([T]) we took 
AGN distances from 3.6 Mpc (the distance to Centaurus 
A) up to 100 Mpc. In the bottom panel of same Fig. ([T]) 
we have plotted the corresponding result for GRBs, at 
typical distances of 10 - 10^ Mpc. The vertical lines cor- 
respond to the relevant region for m2 and ms for the case 
mi = 0. As we will discuss below, putting into a theory 
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context, such couplings are relatively large. In order to f^^^^ ^ using the ratio R,^.,, = ^ 
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where ^12, ^23 and ^13 are the neutrino mixing angles 
determined in neutrino oscihation experiments. Fig. ([2| 
show the expected values for this ratio when the neutrino 
mixing angles lie within the la bands from their current 
global best fit values [26H28] . for which Rv^-.v^ — 8. One 
can also see, Fig. (|3| that by allowing these parameters 
to vary up to their three sigma ranges, Ry^-.v^ can be as 
large as 25 or as low as 2. 
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TABLE I: Ru^-.u^ for current neutrino oscillation fits. 



FIG. 1: Gorrelation between the decaying neutrino mass rrii 
and the required coupling constant gij for the case of AGN's 
(top panel) and GRBs (bottom panel). 



have an estimate of the neutrino flux reduction resulting 
from neutrino decay we note that, since coherence is lost, 
the final flux of a given neutrino flavour will be 
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FIG. 2: Rue-.u^ versus the CP phase 5 for neutrino mixing 
angles at la, with ^23 in the 1** octant (left panel) and 2^^ 
octant (right panel) ]26j . 



where (/)'s are neutrino fluxes at production and detec- 
tion, L the travel distance, the neutrino lifetime in 
the laboratory frame and Uai the elements of the lepton 
mixing matrix [25 . Typical neutrino energies lie in the 
range of 10^ TeV and 10^ TeV for AGNs and GRBs re- 
spectively. Note that, in the limit that L ^ Ti where 
only the stable state survives Eq. ([3| becomes 
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Very similar results are found for global fit of Ref. 




Here we take a normal hierarchy neutrino mass spectrum, 
the disappearance of all states except the lightest (in this 
case vi) is allowed. The final flux of and can 

be computed from Eq. Q and will depend on the three 
mixing angles and the Dirac CP phase 5. In particular, 
we can calculate the suppression of the muon neutrino 



FIG. 3: Rue-.u^ versus the CP phase 6 when the neutrino 
mixing angles lie within 3a of their current best fit value ^6] . 



We now turn to the issue of theoretical consistency of 
the decay hypothesis. In most SU(3)c SU(2)l (8) U(1)y 
seesaw models with spontaneous lepton number violation 
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when one diagonahzes the neutrino mass matrix one also 
diagonalizes, to first approximation, the coupling of the 
resulting Nambu-Goldstone boson to the mass eigenstate 
neutrinos [3 . The exact form of the light-neutrino ma- 
joron couplings can be determined explicitly by pertur- 
bative diagonahzation of the seesaw mass matrix, or by 
using a more general approach using only the symmetry 
properties. The result is [3^ 
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where the subscript S denotes symmetrization, D and 
M are the Dirac and Majorana mass terms in, say, the 
type-I seesaw and Vi is the light neutrino diagonahzation 
matrix. One sees that the majoron couples proportinally 
to the light neutrino mass, hence diagonal to first ap- 
proximation. The off-diagonal part of the gij is inversely 
proportional to three powers of lepton number violation 
scale vi = (cr), since M (x Vi. This is tiny, the only 
hope being to use a seesaw scheme that allows for a very 
low lepton number violation scale, such as the inverse 
seesaw [29]. The particle content is the same as that of 
the Standard Model (SM) except for the addition of a 
pair of two component gauge singlet leptons, and 6'^, 
within each of the three generations, labeled by i. The 
isodoublet neutrinos Ui and the fermion singlets Si have 
the same lepton number, opposite with respect to that 
of the three singlets i/f associated to the "right-handed" 
neutrinos. In the S basis the 9x9 neutral lepton 

mass matrix M., has the form: 
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where rrijj (x (<I>) is the standard Dirac term coming from 
the SM Higgs vev and M is a bare mass term. The term 
/i (X (cr), the vacuum expectation value of a responsi- 
ble for spontaneous low-scale lepton number violation as 
proposed in [7 . This gives rise to a majoron J, 



J = a/2 Imcr. 



(8) 



As a result of diagonahzation one obtains an effec- 
tive light neutrino mass matrix. Note that lepton num- 
ber symmetry is recovered as /i ^ 0, making the three 
light neutrinos strictly massless. The majoron couplings 



of the light mass eigenstate neutrinos are determined 
again as a sum of two pieces as in Eq. (|6|. Detailed 
calculation shows that its off-diagonal part behaves as 
~ Ii^D'^M~^. Even if the M can be significantly 
lower than that of the standard high-scale type-I seesaw 
it is clear that this is way too small in order to produce 
neutrino decay within the relevant astrophysical scales. 

The only way out is to induce a mismatch between the 
neutrino mass basis and the coupling basis. This can be 
achieved by making lepton number a family-dependent 
symmetry [5, 6 . The model is by no means unique, here 
we give an example based on SU(3)c SU(2)l (8) U(1)y 
®U(1)h assigned as 
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TABLE II: Model field content and transformation properties 

) U(1)y ^U{1)h invariant La- 



The SU(3)c SU(2)l 
grangian would be 



yM^..^R^SjCF -^Y^^^SiSj ficTijSiSja, 
where the relevant sub-matrices are 
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One can check explicitly that the first term in Eq. (|6| 
is already non-diagonal and, for sufficiently low values of 
the U{1)h breaking scale can induce a decay sufficiently 
fast as to suppress the flux of to account for the non 
observation of by Ice Cube. 

Finally we comment on a possible indirect test of the 
neutrino decay hypothesis through the observation of 
majoron-emitting JOu^P decay. Indeed if the majoron 
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FIG. 4: Correlating Ru^-.u^ to majoron-emitting JOz//3/3 cou- 
pling {gee)- 

exists and its coupling to the electron neutrino is not 
expected to significantly differ from the one required to 
explain the muon neutrino deficit in IceCube through 
the neutrino decay hypothesis. This correlation is 
depicted in Fig. [4j The figure corresponds to fixing the 
neutrino mixing angles at their best fit values with the 



Dirac mass entry at ttid -^10 GeV, M-1 TeV -1 
TeV, (1 1 KeV and /ia 1 KeV where the ~ sign takes 
into account order one differences among the various 
flavour components of each block. The solid (red) curve 
corresponds to the second octant of the atmospheric 
mixing angle ^23 while the solid (blue) corresponds to 
the first octant, varying the CP phase from to 27r. One 
sees that larged values of the R ratio are preferred by 
the second octant. In conclusion one sees that the decay 
hypothesis invoked to account for the IceCube results 
may be tested in the upcoming searches for the JOiy^/S 
decay. 
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